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ABSTRACT 

We present a set of structural parameters for the central parts of 57 early- 
type galaxies observed with the Planetary Camera of the Hubble Space Telescope. 
These parameters are based on a new empirical law that successfully characterizes 
the centers of early type galaxies. This empirical law assumes that the surface 
brightness profile is a combination of two power laws with different slopes 7 and 
(3 for the inner and outer regions. Conventional structural parameters such as 
core radius and central surface brightness are replaced by break radius r^, where 
the transition between power-law slopes takes place, and surface brightness fib at 
that radius. An additional parameter a describes the sharpness of the break. The 
structural parameters are derived using a x 2 minimization process applied to the 
mean surface brightness profiles. The resulting model profiles generally give very 
good agreement to the observed profiles out to the radius of ~ 10" imaged by the 
Planetary Camera. Exceptions include galaxies which depart from pure power- 
laws at large radius, those with strong nuclear components, and galaxies partly 
obscured by dust. The uncertainties in the derived parameters are estimated 
using Monte-Carlo simulations which test the stability of solutions in the face 
of photon noise and the effects of the deconvolution process. The covariance of 
the structural parameters is examined by computing contours of constant x 2 i n 
multi-dimensional parameter space. 

Subject headings: galaxies: nuclei - galaxies: elliptical and lenticular, cD 

1. Introduction 

Recent Hubble Space Telescope (HST) observations of early-type galaxies reveal that the 
central parts have surface-brightness distributions that are different from the extrapolation 
of traditional fitting formulae derived from ground-based observations (Crane et al. 1993; 
Kormendy et al. 1995; van den Bosch et al. 1994; Forbes, Franx, & Illingworth 1995, 
Lauer et al. 1995, hereafter Paper I). The surface brightness profiles generally consist of 
two distinct regions; a steep power-law regime (/(V) oc r _/3 ) at large radius, and a shallower 
power law (/(r) oc r~ 7 ) at small radius. Galaxies with a quite shallow inner power-law slope 
(7 < 0.3) are classified as "core" galaxies, while those with 7 > 0.5 are labeled power-law 
galaxies (Paper I). Note that 7 is significantly different from zero in most objects studied 
with HST thus far, with only rare exceptions (van den Bosch et al. 1994; Paper I; Ajhar et 
al. 1995). 
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Conventional parameters characterizing the luminosity structure of the central regions 
are the central surface brightness and the core radius, defined as the radius where the surface 
brightness is half the central value. These parameters have been very useful in understand- 
ing elliptical galaxies and bulges of later type galaxies, especially through the fundamental 
plane correlations (Kormendy 1982, 1984, 1985, 1987a, b; Lauer 1985). However, as we do 
not find constant surface brightness cores in most HST observations, we set about finding 
an alternative parameterization (presented below) to better describe the observed surface 
brightness distributions. For a more detailed history on this subject, readers are referred to 
Paper I. 

An alternative approach to the one taken here is that of Gebhardt et al. 1996), who carry 
out a non-parametric analysis of essentially the same sample considered here. The Gebhardt 
et al. analysis is both more thorough and more accurate, but it largely verifies the simple 
fits we describe below. The major advantage of the approach taken in the present paper is 
that it reduces each galaxy to a set of 5 numbers. On the other hand, the non-parametric 
approach does not force galaxies into possibly ill-fitting boxes. 

A one-parameter family of models for spherical stellar systems which yields double 
power-law surface brightness profiles was introduced independently by Dehnen (1993; 7 
models) and by Tremaine et al. (1994; 77 models). In the notation of Tremaine et al. , the 
density 

However, the majority of the actual galaxy profiles is not well-fit by ?7-models. Although 
these models allow for a range of inner power-law slopes, the variations apparent in the 
sharpness of the transition from inner to outer power-laws cannot be reproduced, and the 
fixed asymptotic outer slope of r" 3 does not fit all galaxies. 

To parametrize the HST-observed surface brightness profiles, we introduced in Paper I 
a more general empirical double power-law (the "Nuker" law): 



I(r) = h 2^ (L 



i+i- 

n 



0-1 



(2) 



The break radius, r b , is the radius at which the steep outer profile, J(r) oc r _/3 , "breaks" to 
become the inner shallow profile, I(r) oc r~ 7 . More precisely, is the radius at which the 
absolute logarithmic curvature in the surface brightness profile, \d 2 logI(r)/d(logr) 2 \, is a 
maximum. is the vertical scaling parameter and is simply the surface brightness at r&. As 
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shown in Paper I, the parameter a is necessary to account for the variations in the sharpness 
of the break. 

The Nuker law contains many simpler fitting formulae as special cases; 

• The Hubble-Reynolds law (Reynolds 1913; Hubble 1930), 

(r + a) 1 

corresponds to a — 1, f3 — 2, 7 = 0. 

• The modified Hubble or analytic King law (Binney & Tremaine 1987), 

v c 2 



J(r) = 



r 2 _|_ r 



2 ' 



corresponds to a = 2, (3 = 2, 7 = 0. 

de Vaucouleurs' r 1 / 4 law (de Vaucouleurs 1948), 

I(r) = J exp(-A;r 1/m ), m = 4, 

and its generalization to arbitrary m (Sersic 1968), corresponds to the parameters 
a = ^, 7 = 0, in the limit r fe — > 00, = ^kr\j m — > 00. 

Ferrarese et al. (1994) fit HST photometry of early-type galaxies to the law, 



, r \ —Pi 
I(r) = 2^I C ' ' 



r, 



r N 2(/3 2 -/3i) 

1 + ' 



-1/2 



which corresponds to cc = 2(/3 2 — /3 = fa, 7 = A- 

Galaxies in which the three-dimensional luminosity density varies smoothly near the 
center are said to have "analytic cores" . An analytic core requires that the surface brightness 
near the center has the form I(r) = Iq + I\r 2 + 0(r 4 ). Thus only galaxies with a = 2, 7 = 
have analytic cores. 

The Nuker law also provides excellent fits to the family of 77-models (Figure 1). The 
behavior of the five Nuker parameters for each of the 77-models is illustrated in Figure 2 in 
order to show the parameter space occupied by ^-models. As will be seen below, real galaxies 
cover a much wider range in these parameters. 
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The Nuker law is convenient not only for describing the observed surface brightness 
distribution but also for computing quantities such as the luminosity density and the local 
and line-of-sight velocity dispersion distributions. A FORTRAN code that carries out these 
tasks is available on request from S. Tremaine. 

In this paper, we fit the Nuker law and derive the associated structural parameters for 
57 early-type galaxies observed with the HST. With the exception of NGC 4881, all objects 
were observed prior to the HST refurbishment with the Planetary Camera (0.//043 per pixel) 
using filter F555W (corresponding roughly to Johnson V). The sample includes early-type 
galaxies observed as part of our Guest Observer (GO) programs (#2600, #3912, PI Faber), 
as part of several WFPC Instrument Definition Team programs (#1105, #3229, #3286, 
#3639, & #5233, PI Westphal), and includes F555W PC images of early- type galaxies 
publically available in the HST archive by June 23, 1994. The latter were taken as part 
of GO programs #1038 (PI Ford) and #2607 (PI Jaffe). While small amounts of dust 
could generally be avoided in determining surface brightness profiles, galaxies with severe 
nuclear dust obscuration were rejected. Two or more exposures were used for each galaxy 
to permit cosmic ray rejection, and the peak counts per pixel in the coadded frames were 
typically ~2000 before deconvolution. Further details of data acquisition and reduction of 
the GO #2600/3912 sample are given in Paper I. These reduction procedures have been 
uniformly applied to the additional galaxies in the present sample. The only exception is 
NGC 4881, which was observed after the refurbishment of HST and therefore did not require 
extensive deconvolution. The parameters determined in this paper are used in the analysis 
and interpretation of core properties presented in Faber et al. (1996, Paper III). 



2. Fitting Scheme 

The five structural parameters are determined using a ^-minimization method, x 2 is 
defined as the sum of the squares of the differences in surface brightness, A/i(r) 2 , between the 
observed profiles and those generated by the Nuker law, where /i(r) is in units of magnitudes 
per arcsecond squared. We make use of the mean profiles rather than the major-axis profiles 
published in Paper I. Each galaxy is fitted with a serious of concentric ellipses and the mean 
profile is defined so that for each isophotal ellipse, r = y/r ma j ■ r min , where r ma j and r min are 
major and minor axis radii, respectively. The use of mean profiles rather than major-axis 
profiles reduces the complexities caused by the variation of ellipticity with radius. Significant 
variations in ellipticity may cause the apparent shape of the major axis and mean surface 
brightness distributions to differ quite substantially. 

As discussed in Paper I, the surface brightness measurements were made using two 
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procedures, one optimized for high resolution in the inner arcsecond (Lauer 1985) and the 
other optimized to handle the low signal-to-noise ratio data at larger r. The two profiles were 
subsequently combined. However, the resulting profile has inappropriate spacing between 
data points when plotted in logarithmic space; this is especially true for the profile obtained 
from the high resolution optimized technique, for which the point spacing is linear. In order 
to treat all radii with equal importance for our power-law fitting, the observed profiles were 
resampled beginning at the first measured data point (0.044"/ \^2n = 0.017") and progressing 
at equal logarithmic intervals of Alogr = 0.05 using spline interpolation. This spacing is 
small enough to represent the variation of surface brightness in the original profiles. Uniform 
weights and errors were applied to the individual data points in the resampled profile during 
the fitting procedure. 

For the minimization we have used the MINUIT package, which handles non-linear 
minimization through a combination of Monte-Carlo search, downhill-Simplex, and the vari- 
able metric method (James and Roos 1977). Due to the optical problems inherent in pre- 
refurbishment HST data and the high central concentration of light in the galaxies them- 
selves, we applied the following additional adjustments to the minimization procedure. 

Even with an accurate point-spread function (PSF) and careful image deconvolution, 
the effect of the spherical aberration of the HST primary cannot be completely corrected 
within 0."1 of the galaxy centers (Lauer et al. 1992b). Therefore the fitting domain is 
restricted to radii greater than 0."1. Nevertheless profiles do differ greatly within 0."1, so 
there is some information there. To utilize this, the total light within 0."1 (as opposed to its 
distribution) was used as an additional constraint on the the fit. The weight of the central 
light is equivalent to the sum of weights for the resampled data points which lie within 0."1. 
Note that the constraint given by the central light does not necessarily force the fit to pass 
through the data points within 0."1. Large departures between the fit and the innermost 
data point often occur. This is because the total light is easily adjusted by only a slight 
change in the surface brightness in the outer region (but still within 0."1), where the areal 
coverage is comparatively large. 

Nuker law fits to the sample galaxies are shown in Figure 3, and the structural param- 
eters are listed in Table 1. For a few galaxies with nuclei such as NGC 1331, NGC 3599, 
and NGC 4239, a separate fit was made to the bulge component only, excluding the nucleus. 
In these cases, the central light constraint is also ignored in the fit. Any dust patches were 
masked out from the raw image before the profile extraction. Nevertheless the galaxies with 
significant dust features, such as NGC 524, show distinct distortions in their profiles. In 
these cases, we set the weights to zero in the regions obviously affected by dust obscuration. 
For each galaxy, the fitting domain is shown by the solid symbols. The RMS deviation 
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computed over the fitting domain between the model and the original surface brightness 
measurements is listed in Table 1. 

While the Nuker fits generally match quite well the observed profiles out to the outermost 
HST data points (r ~ 10"), there are cases like Abell 1020, NGC 4478, and NGC 4486B, 
where the outer parts deviate (both upward and downward) from the Nuker law. Even 
galaxies which show good agreement with the Nuker law within 10" in general will also fail 
at much larger radii beyond the field covered by the present HST data, as the profiles follow 
a curving de Vaucouleurs law, not a power law there. For some of the sample galaxies, 
the HST data are combined with ground based photometry covering much larger area and 
presented in Kormendy et al. (1996, in preparation). 

3. Parameter Uncertainties 

3.1. Error Estimates 

The deconvolution process can introduce systematic errors into surface-brightness pro- 
files, particularly if use is made of an inappropriate PSF. Experiments with different PSFs 
were carried out in Paper I and revealed both photometric offsets and changes in slope out to 
~ 2". These effects are more pronounced for galaxies with steep surface brightness profiles 
near the center. However, the experiments showed that even for galaxies with particularly 
sharp nuclei, the systematic errors in the inner parts (where the effects of deconvolution are 
greatest) are small (~ 0. m 06). We note that the slope change would be echoed primarily by 
7, a, and h; (3 and would be essentially unaffected. 

The deconvolution process also amplifies photometric uncertainties, and we have carried 
out simulations to better understand the errors in the Nuker parameters introduced by 
deconvolution of data with a finite signal-to-noise ratio. Each simulation was carried out 
as follows: (i) Poisson noise appropriate to the number of collected photons in each pixel 
was added to the raw (undeconvolved) image of a galaxy, (ii) Each such realization was 
deconvolved using 80 iterations of the Lucy- Richardson deconvolution algorithm. (Hi) The 
two-dimensional data were fitted with a series of concentric ellipses as described in Paper I. 
(iv) The mean profiles were fitted with the Nuker law and the resulting parameters stored. 

It is important to note that we have added artificial noise on top of the natural noise 
already present in the data. In this sense, our simulations overestimate (by y/2) the errors 
due purely to S/N. To isolate the effects of S/N alone, we generated a noiseless model image 
and repeated the experiments. By comparing the results of these experiments to those using 
real data, we were able to estimate to what extent the goodness of fit was determined by 
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non-photon statistical errors, such as changes in ellipticity, flat-field defects, unsubtracted 
stars, undetected dust obscuration, etc. 

We chose three galaxies for this exercise : Abell 1831, NGC 4889, and NGC 4434. Abell 
1831 and NGC 4889 are typical examples of well resolved core galaxies with small 7, while 
NGC 4434 is a power-law galaxy. These galaxies also span a large range in central signal- 
to-noise ratio. Abell 1831 has the lowest central signal-to-noise ratio in the sample, with a 
peak data count around 1.0 x 10 3 DN (digital numbers) per pixel. By comparison, NGC 4889 
has a peak count of 6.9 x 10 3 DN per pixel, and NGC 4434 has a peak count of 1.1 x 10 4 DN 
per pixel. 

For each galaxy one hundred realizations of the original image were generated. These 
realizations were then deconvolved and independently analyzed. We constructed a symmet- 
ric, two-dimensional, noiseless model image based on the best-fit parameters for Abell 1831, 
convolved with an appropriate PSF, and scaled to match the signal-to-noise ratio in the real 
image. Experiments with this model image provided a worst-case upper limit on the degree 
to which the results in Table 1 are affected by purely statistical uncertainties. 

The profiles fitted to the deconvolved, simulated images are shown in Figure 4. As 
expected, Abell 1831 and its model representation show the largest spread among the three 
galaxies, and deviations from the best-fit profile are apparent out to r ~ 1". NGC 4889 
also shows an obvious dispersion, but in this case the deviations are generally confined to 
r < 0."1, and the fit parameters are rarely significantly different from their best-fit values. 
NGC 4434 has the highest signal-to-noise ratio and, despite the steepness of the profile near 
the center, the effects of photon noise and deconvolution appear to be minimal. 

We note that deconvolving a natural image with added noise changes the shape of the 
simulated profiles somewhat, at least in the case of Abell 1831 and its model represenation. 
There is a slight increase in surface brightness in the region 0."1 < r < 1". The Nuker 
parameters derived from the simulated profiles therefore are somewhat different from those 
listed in Table 1. This effect is not detected in the other two galaxies with higher signal- 
to-noise ratios, and we attribute it to the non-conservative nature of the Lucy-Richardson 
deconvolution algorithm employed. 

The frequency distributions of the best-fit Nuker parameters for the four galaxy images 
are plotted in Figure 5. The mean parameter value and associated standard deviation in 
each panel were determined by fitting a Gaussian to each frequency histogram. As expected 
from the amplitude of the variations apparent in the simulated profiles, Abell 1831 shows 
the largest standard deviation for each of its Nuker parameters except (3. (3, the slope of the 
outer power-law, is rather stable in Abell 1831 because r b is comparatively small and the 
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outer power law portion of the profile is sampled over a large radial range. 

The average values of the Nuker parameters for Abell 1831 and its model representation 
are also not identical. This is due in part to the different noise characteristics of the two 
images near the break radius and their amplification in the deconvolution process. It is also 
due to the fact that the model image does not exactly mimic the real data. For example, 
the slight dip apparent at r ~ 0."15 in the real data is not reproduced in the model image, 
and the best fits will differ accordingly. 

We note that there is no significant difference in x 2 between the best-fit Nuker pa- 
rameters for Abell 1831 and those for its noiseless model representation. This implies that 
our computed x 2 is much less dependent on signal-to-noise ratio than it is on other factors. 
These factors would include image defects as well as real discrepancies in form between Nuker 
profiles and those of galaxies. 

In spite of its low signal-to-noise ratio, both Abell 1831 and its model representation 
show that errors in the Nuker parameters introduced by photon noise and deconvolution 
are small. We conclude that the entries in Table 1 are largely unaffected by differences in 
signal-to- noise ratio. 

3.2. Covariances 

The uncertainties in the derived parameters and the coupling of the parameters with 
one another are usually investigated using the distribution of x 2 ■ For any two parameters, 
we can project the 5-dimensional \ 2 distribution onto the appropriate axes. This reveals 
not only the how individual parameters are constrained by the data, but also the degree to 
which various parameters are correlated. However, since the data have been resampled and 
assigned uniform weights in order to put equal emphasis on different portions of the surface 
brightness profiles, the values of \ 2 we compute, while relatively correct, are not properly 
normalized to compute absolute uncertainties. Our use of the total light within 0."1 as an 
arbitrarily weighted, additional constraint exacerbates this problem. We could conceivably 
compute an a priori x 2 based on photon statistics, but that would not include errors due to 
flat-fielding defects, improper deconvolution, and systematic mismatch betweeen the Nuker 
law and the galaxies. 

In Figure 6 we show projected, constant-^ 2 contours for several combinations of Nuker 
parameters for the three galaxies considered in the previous section. For the reasons given 
above, the confidence limits for each parameter cannot be directly estimated from the initially 
computed x 2 values. To remedy this, we have normalized the computed x 2 values using the 
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observed distributions of best-fit parameters found in the Monte Carlo simulations. The 
inner, thick contours in Figure 6 thus correspond to the 1 a errors determined from the 
Monte Carlo simulations. In order to show the gradient of the projected x 2 distribution, the 
outer, thick contours are used to indicate the 10 a level. This Monte Carlo normalization 
includes photon statistics but not the other error sources note above, so the treatment is 
approximate. 

The general conclusion from these plots is that, while certain error pairs are correlated, 
the only truly serious correlation is between and /x& for power law galaxies. This is to be 
expected, as these parameters are formally completely unconstrained for pure power laws 
with zero curvature. 

If we exclude those galaxies with nuclei (NGC 1331, 3599, 4239, and 6166) or with 
obvious and widespread dust (NGC 524), the remaining 52 galaxies give a mean RMS de- 
viation between the best-fit Nuker law and the surface brightness measurements of 2.7% 
(0.03 mag). This is comparable to, though somewhat larger than, an RMS of 1% obtained 
by Gebhardt et al. (1996) using smoothed spline fits to the data. This simply reflects the 
fact that the Nuker law is a smoother estimate of the surface brightness compared with the 
non-parametric fit. For the purposes of determining the physical properties of an individual 
galaxy, it is therefore preferable to invert a non-parametric form of the surface brightness 
profile. However, if one is searching for global trends and systematics, the Nuker param- 
eters provide a convenient basis for comparison. The conclusions reached in Paper I and 
Faber et al. (1996) based on the Nuker parameters are not altered with the application of 
non-parametric methods. 



4. Summary 

In this paper, we derive a set of structural parameters for 57 early-type galaxies observed 
with HST. These parameters are based on the empirical "Nuker" law, which successfully 
describes the central surface-brightness distributions within the radius range from ~0."1 to 
~10" for most galaxies in the sample. These parameters have been used in Paper I and Faber 
et al. (1996) to identify distinct classes of objects, characterized by significant differences in 
one or more of the Nuker parameters. While these parameters can in principle also be used to 
evaluate luminosity densities and (with the availability of suitable spectral information) the 
run of velocity dispersion with radius, we believe this would be better left to non-parametric 
methods (Gebhardt et al. 1996), which are more faithful to small but significant undulations 
in the observed profiles. 
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Fig. 1. — Surface brightness profiles for ?7-models with a range of i] values (open circles). 
Solid line represents Nuker law fits to each model profile. 

Fig. 2. — Nuker parameters for 77 models. Left diagrams are for the fits made to the radius 
range from 1CT 5 to 10 3 . Right diagrams are with the restricted radius range of 0.1 < r < 10, 
which is comparable to our fitting regime for the sample galaxies, fib is the surface brightness 
at r b , i.e. fi b = -2.51og 10 i&. 

Fig. 3. — Mean surface brightness profiles (filled and open circles) and Nuker law fits (solid 
lines) for the sample galaxies. The surface brightness is in mag /sec 2 . The filled circles 
indicate the domain over which the fits are were computed, while the open circles denote 
the points excluded from the fit. The data points inside r ~ 0."1 are regarded as uncertain, 
therefore are not directly used in the fitting (see text). 

Fig. 4.— Simulated surface brightness profiles (solid lines) for Abell 1831, NGC 4889 
and NGC 4434, shown with the original profile (open circles). The top diagram is for the 
noiseless model image of Abell 1831. The peak data counts are 1.0 x 10 3 DN for Abell 1831, 
6.9 x 10 3 DN for NGC 4889, and 1.1 x 10 4 DN for NGC 4434 respectively. The means and 
RMSs of the numerical experiments at 0."03 and 0."05 are shown as open squares with error 
bars. 

Fig. 5. — Frequency histogram of Nuker parameters derived from 100 simulated profiles 
of the noiseless model image of Abell 1831, Abell 1831, NGC 4889, and NGC 4434. A 
Gaussian distribution is evaluated for each parameter and shown by solid curve. The mean 
and standard deviation values are given, fib is in mag /sec 2 and is in arcsec. 

Fig. 6. — Parameter cross-correlations shown by the x 2 ellipses for Abell 1831, NGC 4434, 
and NGC 4889. fi b is in mag /sec 2 and r b is in arcsec. The thick solid curves indicate la 
and 10a ellipses as scaled by the results of Monte-Carlo simulations given in Figure 5. The 
locations of the minimum x 2 are indicated by crosses. 




Radius 



TABLE 1 : STRUCTURAL PARAMETERS 



Galaxy 


Sample 


4 


n 


a f 


3 7 


RMS (%) 


Note 



On 
On 

PL, 

(N 
(N 



> 



(N 

o 

Oh: 
6 



ADCll 1UZU 


oylZ 


171 an 
1 1 .luU 


U.Zou 


Z.0D4 


1 QO/1 

i.oy4 


n 1 7E. 
U. 1 ( o 


1 Q/l 

l.o4 




1 Q Q 1 "f 

ADCll looli 


2600 


18.698 


0.500 


3.567 


1.170 


0.110 


2.19 




Abcll ZUbz 


2600 


18.600 


0.419 


8.022 


0.749 


0.201 


3.71 


• 

nuclear point source excluded in fit 


IN LxO oz4 


3912 


16.121 


0.322 


1.293 


1.001 


0.005 


5.36 


dust region excluded in fit 


tvt/^^i cnft 


2600 


18.116 


3.501 


0.758 


1.968 


0.554 


3.37 




IN CrO f ZU 


2600 


17.500 


3.208 


2.315 


1.660 


0.062 


0.51 




lNi_xL< lUZo 


2600 


16.357 


1.865 


4.722 


1. 184 


0.782 


2.28 




Mfi p ii 70 


3912 


18.694 


2.431 


1.517 


1.638 


1.009 


3.43 




f" 1 1 QQ 1 

1N<jU lool 


3912 


21.035 


8.157 


1.378 


2.842 


0.472 


2.22 


nucleus excluded in fit ( bulge ) 


1N<jU lool 


3912 


19.976 


4.083 


4.470 


1.615 


0.673 


5.76 


nucleus included in fit 


iN<jU loyy 


2600 


17.059 


3. 144 


1.504 


1.682 


0.068 


0.70 




INOU 14UU 


3912 


15.510 


0.330 


1.395 


1.319 





3.20 




AT /~l /~1 1 /I Oft 

INOU 14zo 


3912 


17.540 


1.644 


3.619 


1.354 


0.847 


3.06 




JNCjC loUU 


3639 


18.951 


5.133 


1.255 


2.177 





2.27 




Mrin i '7(\r\ 
INCjL/ 1/UU 


3912 


14.041 


0.088 


0.903 


1.303 





2.39 




Oft Q ft 

INCjO ZooO 


3912 


15.708 


0.093 


1.838 


1. 140 


0.042 


3.11 




at t~< o Q Q O 

INLxO ZooZ 


3912 


17.452 


0.910 


1.837 


1.399 


0.024 


1.05 




INLxO Zo41 


3912 


14.548 


0.133 


0.928 


1.021 


0.009 


3.55 




MPf oil c 

INLxL/ olio 


oftnn 
ZOUU 


1 d OA Q 

lD.Z4c5 


z.yuo 


1 A ftfl 

1.4DD 


i a on 

i.4zy 


n 701 

u. < yi 


i en 

i.oy 




Mrin QQ77 

inCjO oo/ r 


2600 


12.902 


0.089 


1.923 


1.332 


0.292 


2.71 




AT (~1 QQ TO 

INCjO oo/y 


3229 


16.136 


1.738 


1.593 


1.432 


0.182 


2.05 




in^jLj ooyy 


oylZ 


1 1 Qnn 
1 1 .oyu 


i oi n 

i.ziy 


i q 7ft n 
lo. f DU 


i ft7n 
l.D ( U 


n con 
U.ooU 


7 A Q 

/ .4o 


nucleus excluded in fit ( bulge ) 


in^xLj ooyy 


QOI O 

oylZ 


1 1 t;7Q 
1 / .o / o 


1 Q C 1 

l.ool 


i q nnn 

lo.uuy 


i ftftn 
l.DDU 


n 70/i 
U. ( y4 


n i o 
y. IZ 


nucleus included in fit 


Mnn Qftnc 
INLxLj oOUo 


oftnn 
ZOUU 


1 1 QK.Q 
1 I .ZOO 


U.ooD 


y.ioo 


1 OftQ 

l.ZDo 


n ft7i 
U.D ( 1 


O it; 
Z. lo 




Mfp QftOQ 

INLxL^ oOUo 


oylZ 


1 c; A A O 

io.44y 


n oq /i 
U.ZS4 


l.UoZ 


i QQn 
l.ooU 


n 
U 


1 C E" 

l.oo 




IN uv 4:100 


ZDU / 


1 Q Q K7 
lo.OO / 


9 ^91 
Z. OZl 


u.y^o 


i ^nn 

1 . OUU 


n 1 ac\ 


1 ^ 
1.00 




IN 4zoy 


QOI 9 

oy iz 


1 Q f!9S 
lo.UZo 


U.olo 


A OQ7 
4.ZO / 


n nft7 
u.yD ( 


n qe;o 
U.ooZ 


q on 

O.ZU 


nucleus excluded in fit ( bulge J 


IN 4zoy 


QOI O 

oy IZ 


lo.4ZU 


1 9Q/1 
1 .Zo4 


1 A f^QI 
14.0O1 


n n^n 
u.yoy 


n ft/ift 

U.D4D 


7 fll 

/ .ol 


nucleus included in fit 


In 4Z0o 


lUoo 


1 Q 71 9 
lO. / IZ 


U.U40 


1 fl/1 7 
1.U4 ( 


n Q7Q 

U.o/o 


n nno 

U.UUZ 


Q 7/1 
O. ( 4 




tvt CT* /IQftr; 
INt^xO 4oOO 


ZDU / 


1 K A A C 
10.44O 


1 A71 
1.0(1 


Z.UoD 


i 07n 
l.Z / u 


n 1 /i a 

U. 14o 


"1 79 
1. ( Z 




IN 4oo / 


oftnn 
ZoUU 


1 Q OOQ 

lo.yyo 


-1 /I OQ 
4.4Zo 


Q QftO 

o.oDZ 


i c;m 

i.oyi 


n 7i n 

u. < iy 


O C 1 
Z.ol 




IN 44o4-|- 


oftnn 
ZDUU 


1 7 Q /I O 

1 / .o4Z 


1 OQD 

l.ZoU 


o ono 

z.suy 


1 QQ c; 
l.ooo 


n 7(1(1 

u. < yy 


O Qft 

Z.oD 




IN LxL/ 44oo 


oftnn 
ZDUU 


1/1 cr /I "I 
14.041 


n no 
U. llo 


c: oft c; 
o.ZDo 


1 A 07 
1.4Z ( 


n /i nn 

u.4yu 


O QQ 

Z.oo 




IN LxL/ 4404 


oftnn 
ZDUU 


1 f .o4o 


1 OI o 
1.Z1Z 


1 ft/l o 
1.D4Z 


1 ftQQ 

l.Doo 


n QQ A 
U.oo4 


i nt; 

i.yo 




IN LxL/ 440 ( 


oftnn 
ZDUU 


i o noo 

ly.yoz 


Q OA Q 

o.Z4o 


7 El 

/ .olo 


O 1 QQ 

Z. loo 


n noi 

u.ysi 


c; ftQ 
o.Do 




Mfp AAVO 

IN LxL/ 44 / Z 


QOOO 

ozzy 


1 ftc;o 

lo.Doy 


O QOQ 

Z.oyo 


o nQi 
Z.Uol 


1 1 ftC 


n nQn 

u.uoy 


1 Q1 

l.ol 




\Tpp A ATI 

IN LxL/ 44 / o 


oftn7 
ZDU / 


17 QQC 

1 f .000 


Q "7C\A 

o. /U4 


n ftft k 
U.DDo 


O A A 7 
Z.44 ( 


n m a 
U.U14 


i no 

i.yz 


double nucleus : 


MP A A 7Q 
IN LxLj 44 ( o 


oftn7 
ZDU f 


10.400 


n 1 7K. 

u. i ( o 


Q QOI 

o.oZl 


n QQft 
U.ooD 


n /io7 
U.4Z / 


O C A 
Z.o4 




wpp /1/1Qft 
IN VjU 44o0 


i in^ 
llUo 


1 7 OOQ 

i ( .yzo 


( .DUD 


O S1 ft 
Z.olD 


1 QO/1 

i.oy4 


n o/ift 
U.Z4D 


n fto 

U.DZ 


point source removed; profile from Lauer Gt 0,1 (1992a) 


A A Qft"D 


oyiz 


1 /I 07ft 

14. y ( d 


n i oi 
U. lol 


O 7Q1 

Z. / ol 


i Qon 

i.ozy 


n ion 

u. ioy 


1 OQ 
l.Zo 




INijxO 400 1 


oftnn 
ZDUU 


1 Q OOO 

lo.yzz 


Q QftO 

O.oDZ 


O CiA A 
Z.y44 


1 OQQ 

l.Zoo 


n ono 
U.oUZ 


o ftt; 
Z.Do 




INijxO 4ooZ 


QOQft 

oZoO 


1C KIO 

lo.oiz 


n a a n 
U.D4y 


1 /ICQ 

1.4oo 


1 007 

i.zy ( 


n 
U 


A Q7 
4.0 ( 




Mpp /ICft/l 
INijxO 4004 


ZDU / 


1 c; 700 
lO. ( ZZ 


U.olD 


n oc; c; 
U.ZDO 


i nnc; 

i.yuo 


n n/i 7 
U.U4 ( 


o.oo 




IN V_xO 4:0 1 U 


9R07 


1 7 9Q9 


9 Q 
z.oiy 


71 ^ 
O. I -LO 


1 491 


U.O^iO 


2 31 




MfP Aft.01 
1Nv_x^_j 'iOZl 


ozzy 


1 7 9E.fl 


9 Q99 

z.yzz 


n i S7 

U. lo I 


i 7i n 

JL. I -LU 


u.'rtyo 


9 it; 
0.10 




MfP /lftQft. 

1Nv_x^_j 4:000 


oy iz 


1 7 7Q9 
1 1 . i OZ 


^ 91 
0. Z lO 




1 'iOQ 

i .ozy 


n i oo 
u. izy 


1 ^fi 
1 .OO 




NGC 4649 


3229 


17.191 


3.584 


2.005 


1.300 


0.153 


2.62 




NGC 4697 


3912 


16.957 


2.575 


24.863 


1.040 


0.742 


4.12 


dust disk excluded from fit 


NGC 4742 


2600 


16.764 


1.386 


48.599 


1.994 


1.091 


4.64 




NGC 4874 


3912 


19.222 


2.632 


2.327 


1.366 


0.134 


1.23 




NGC 4881 


5233 


17.330 


0.498 


1.763 


1.358 


0.757 


3.25 


post refurbishment data 


NGC 4889J 


3912 


18.053 


1.675 


2.613 


1.348 


0.047 


1.26 




NGC 5813 


3912 


16.529 


0.792 


2.147 


1.326 


0.076 


1.66 




NGC 5845 


2600 


17.625 


2.273 


1.266 


2.742 


0.515 


7.61 




NGC 6166 


3912 


19.348 


2.221 


3.320 


0.992 


0.081 


1.18 


nuclear point source excluded in fit 


NGC 7332 


2600 


15.803 


0.767 


4.250 


1.343 


0.897 


3.24 




VCC 1199 


2600 


19.677 


1.500 


7.986 


1.619 


1.129 


3.99 


VCC galaxies from Binggclli et al (1985) 


VCC 1440 


2600 


19.994 


2.522 


5.537 


1.583 


0.964 


5.17 




VCC 1545 


2600 


19.710 


1.205 


7.652 


1.018 


0.625 


4.69 




VCC 1627 


2600 


20.153 


2.989 


2.123 


2.099 


0.948 


2.37 





f /!/, = — 2.51og q Q/f, ; mag/arcsec 2 

| galaxies used for detailed analysis in this paper 



10" 5 < r < 10 3 



0.1 < r < 10 



4 - 



a. 



3 - 



• • • • • 

J I I I I I I I I I L 



1 .2 

□ 

0.8 
1 .2 

0.9 



J I I I I I I I I I L 



Jl 



J I I I I I I I I I L 



3.1 
^ 3 
2.9 



J I I I I I I I I I L 



^ 1 " 



tt_i i i I i i i I i * i * * 



• • • 



I I I I I I I I I I I 



I I I I I I I I I I I 



• • • 



I I I I I I I I I I I 



• . • . • • . 

• • • 



J I I I I I I I I I L 



J i i i I i i * * T~ 



J I I I I I L 



12 3 

7] 



12 3 

T) 




Mean Radius (arcsec) 



MODEL-Abell 1 831 




18.4 18.6 
AbelM 831 



0.4 0.5 



1.16 1.18 0.1 0.2 



T 



30 
20 
10 




<M„> 18.674 
a 0.02 



± 



± 



<r b > 0.513 
a 0.015 



J_ 



± 



<n> 3.764 
a 0.445 




</?> 1.193 
a 0.007 




<7> 0.109- 
a 0.028 




18.6 18.8 
NGC4889 



0.5 0.6 



1.15 1.2 



0.1 0.2 




18 18.1 
NGC4434 




17.32 17.34 
Mb 



0.75 0.8 
7 



Abell 1 831 




NGC4889 




NGC4434 




